A vortex-based model for the subgrid flux of a passive scalar

A model for the flux of a passive scalar by the subgrid motions in the large-eddy simulation of turbulent flow is proposed within the framework of the stretched-vortex subgrid stress model. The model is based on an analytical solution for the winding of a scalar field by an elemental subgrid vortex. This gives a tensor gradient-diffusion expression for the local flux of the scalar with subgrid turbulent diffusivity which depends upon the subgrid energy, the local cell size, and the vortex orientation in space. The scalar-flux subgrid model is tested by comparison of the results of 323 large-eddy simulation of passive-scalar transport by forced isotropic turbulence in the presence of a mean scalar gradient, with the direct-numerical simulation results of Overholt and Pope [Phys. Fluids 8, 2128 (1996)]. The present large-eddy simulation results predict that at large Taylor–Reynolds numbers, the ratio of the scalar variance to the squared product of the scalar gradient with the dissipation length of the turbulence, is asymptotic to a nearly constant value c[prime]2/(alpha1 Lepsilon)2[approximate]0.36

Pressure spectra for vortex models of fine-scale homogeneous turbulence

Pressure spectra at large wave numbers are calculated for Lundgren–Townsend vortex models of the fine scales of homogeneous turbulence. Specific results are given for the Burgers vortex and also for the Lundgren-strained spiral vortex. For the latter case, it is found that the contribution to the shell-summed spectrum produced by the interaction between the axisymmetric and nonaxisymmetric components of the velocity field is proportional to k^–7/3 (k=||k|| is the modulus of the wave number) in agreement with Kolmogorov-type dimensional arguments. Numerical estimates of the dimensionless prefactors for this component are obtained in Kolmogorov scaling variables and comparisons are made with results from the Batchelor–Kolmogorov theory, and with experimental measurement

Contour Dynamics Methods

In an early paper on the stability of fluid layers with uniform vorticity Rayleigh (1880) remarks: "... In such cases, the velocity curve is composed of portions of straight lines which meet each other at finite angles. This state of things may be supposed to be slightly disturbed by bending the surfaces of transition, and the determination of the subsequent motion depends upon that of the form of these surfaces. For co retains its constant value throughout each layer unchanged in the absence of friction, and by a well-known theorem the whole motion depends upon [omega]." We can now recognize this as essentially the principal of contour dynamics (CD), where [omega] is the uniform vorticity. The theorem referred to is the Biot-Savart law. Nearly a century later Zabusky et al (1979) presented numerical CD calculations of nonlinear vortex patch evolution. Subsequently, owing to its compact form conferring a deceptive simplicity, CD has become a widely used method for the investigation of two-dimensional rotational flow of an incompressible inviscid fluid. The aim of this article is to survey the development, technical details, and vortex-dynamic applications of the CD method in an effort to assess its impact on our understanding of the mechanics of rotational flow in two dimensions at ultrahigh Reynolds numbers. The study of the dynamics of two- and three-dimensional vortex mechanics by computational methods has been an active research area for more than two decades. Quite apart from many practical applications in the aerodynamics of separated flows, the theoretical and numerical study of vortices in incompressible fluids has been stimulated by the idea that turbulent fluid motion may be viewed as comprising ensembles of more or less coherent laminar vortex structures that interact via relatively simple dynamics and by the appeal of the vorticity equation, which does not contain the fluid pressure. Two-dimensional vortex interactions have been perceived as supposedly relevant to the origins of coherent structures observed experimentally in mixing layers, jets, and wakes, and for models of large-scale atmospheric and oceanic turbulence. Interest has often focused on the limit of infinite Reynolds number, where in the absence of boundaries, the inviscid Euler equations are assumed to properly describe the flow dynamics. The numerous surveys of progress in the study of vorticity and the use of numerical methods applied to vortex mechanics include articles by Saffman & Baker (1979) and Saffman (1981) on inviscid vortex interactions and Aref (1983) on two-dimensional flows. Numerical methods have been surveyed by Chorin (1980), and Leonard (1980, 1985). Caflisch (1988) describes work on the mathematical aspects of the subject. Zabusky (1981), Aref (1983), and Melander et al (1987b) discuss various aspects of CD. The review of Dritschel (1989) gives emphasis to numerical issues in CD and to recent computations with contour surgery. This article is confined to a discussion of vortices on a two-dimensional surface. We generally follow Saffman & Baker (1979) in matters of definition. In two dimensions a vortex sheet is a line of discontinuity in velocity while a vortex jump is a line of discontinuity in vorticity. We shall, however, use filament to denote a two-dimensional ribbon of vorticity surrounded by fluid with vorticity of different magnitude (which may be zero), rather than the more usual three-dimensional idea of a vortex tube. The ambiguity is unfortunate but is already in the literature. Additionally, a vortex patch is a finite, singly connected area of uniform vorticity while a vortex strip is an infinite strip of uniform vorticity with finite thickness, or equivalently, an infinite filament. Contour Dynamics will refer to the numerical solution of initial value problems for piecewise constant vorticity distributions by the Lagrangian method of calculating the evolution of the vorticity jumps. Such flows are often related to corresponding solutions of the Euler equations that are steady in some translating or rotating frame of reference. These solutions will be called vortex equilibria, and the numerical technique for computing their shapes based on CD is often referred to as contour statics. The mathematical foundation for the study of vorticity was laid primarily by the well-known investigations of Helmholtz, Kelvin, J. J. Thomson, Love, and others. In our century, efforts to produce numerical simulations of flows governed by the Euler equations have utilized a variety of Eulerian, Lagrangian, and hybrid methods. Among the former are the class of spectral methods that now comprise the prevailing tool for large-scale two- and three-dimensional calculations (see Hussaini & Zang 1987). The Lagrangian methods for two-dimensional flows have been predominantly vortex tracking techniques based on the Helmholtz vorticity laws. The first initial value calculations were those of Rosenhead (193l) and Westwater (1935) who attempted to calculate vortex sheet evolution by the motion of O(10) point vortices. Subsequent efforts by Moore (1974) (see also Moore 1983, 1985) and others to produce more refined computations for vortex sheets have failed for reasons related to the tendency for initially smooth vortex sheet data to produce singularities (Moore 1979). Discrete vortex methods used to study the nonlinear dynamics of vortex patches and layers have included the evolution of assemblies of point vortices by direct summation (e.g. Acton 1976) and the cloud in cell method (Roberts & Christiansen 1972, Christiansen & Zabusky 1973, Aref & Siggia 1980, 1981). For reviews see Leonard (1980) and Aref (1983). These techniques have often been criticized for their lack of accuracy and numerical convergence and because they may be subject to grid scale dispersion. However, many qualitative vortex phenomena observed in nature and in experiments, such as amalgamation events and others still under active investigation (e.g. filamentation) were first simulated numerically with discrete vortices. The contour dynamics approach is attractive because it appears to allow direct access, at least for small times, to the inviscid dynamics for vorticity distributions smoother than those of either point vortices or vortex sheets, while at the same time enabling the mapping of the two-dimensional Euler equations to a one-dimensional Lagrangian form. In Section 2 we discuss the formulation and numerical implementation of contour dynamics for the Euler equations in two dimensions. Section 3 is concerned with applications to isolated and multiple vortex systems and to vortex layers. An attempt is made to relate this work to calculations of the relevant vortex equilibria and to results obtained with other methods. Axisymmetric contour dynamics and the treatment of the multi-layer model of quasigeostrophic flows are described in Section 4 while Section 5 is devoted to a discussion of the tendency shown by vorticity jumps to undergo the strange and subtle phenomenon of filamentation

Vortex ring formation at tube and orifice openings

The formation, at tube and orifice openings, of vortex rings generated by a piston moving with velocity proportional to time to some power m, is considered. The expansion of the axisymmetric generating flow about the circular forming edge is used in conjunction with the similarity theory of edge vortex growth to model the ring formation process. For large Reynolds numbers the ring diameter and circulation are not strongly dependent on the piston velocity profile. However, the ring viscous subcore shows peaks in the tangential velocity profile only if m < (π–θ)/(2π–θ), where θ is the edge forming angle

Direct numerical simulation and large-eddy simulation of stationary buoyancy-driven turbulence

We report direct numerical simulation (DNS) and large-eddy simulation (LES) of statistically stationary buoyancy-driven turbulent mixing of an active scalar. We use an adaptation of the fringe-region technique, which continually supplies the flow with unmixed fluids at two opposite faces of a triply periodic domain in the presence of gravity, effectively maintaining an unstably stratified, but statistically stationary flow. We also develop a new method to solve the governing equations, based on the Helmholtz–Hodge decomposition, that guarantees discrete mass conservation regardless of iteration errors. Whilst some statistics were found to be sensitive to the computational box size, we show, from inner-scaled planar spectra, that the small scales exhibit similarity independent of Reynolds number, density ratio and aspect ratio. We also perform LES of the present flow using the stretched-vortex subgridscale (SGS) model. The utility of an SGS scalar flux closure for passive scalars is demonstrated in the present active-scalar, stably stratified flow setting. The multi-scale character of the stretched-vortex SGS model is shown to enable extension of some second-order statistics to subgrid scales. Comparisons with DNS velocity spectra and velocity-density cospectra show that both the resolved-scale and SGS-extended components of the LES spectra accurately capture important features of the DNS spectra, including small-scale anisotropy and the shape of the viscous roll-off

Large-eddy simulation and wall modelling of turbulent channel flow

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Kármán-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Re_τ based on the friction velocity u_τ and the channel half-width δ in the range 2 × 10^3 to 2 × 10^7. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jiménez (2006) at Re_τ = 2003 and maintain an O(1) grid dependence on Re_τ

A vortex-based subgrid stress model for large-eddy simulation

A class of subgrid stress (SGS) models for large-eddy simulation (LES) is presented based on the idea of structure-based Reynolds-stress closure. The subgrid structure of the turbulence is assumed to consist of stretched vortices whose orientations are determined by the resolved velocity field. An equation which relates the subgrid stress to the structure orientation and the subgrid kinetic energy, together with an assumed Kolmogorov energy spectrum for the subgrid vortices, gives a closed coupling of the SGS model dynamics to the filtered Navier-Stokes equations for the resolved flow quantities. The subgrid energy is calculated directly by use of a local balance between the total dissipation and the sum of the resolved-scale dissipation and production by the resolved scales. Simple one- and two-vortex models are proposed and tested in which the subgrid vortex orientations are either fixed by the local resolved velocity gradients, or rotate in response to the evolution of the gradient field. These models are not of the eddy viscosity type. LES calculations with the present models are described for 32^(3) decaying turbulence and also for forced 32^(3) box turbulence at Taylor Reynolds numbers R-lambda in the range R(lambda)similar or equal to 30 (fully resolved) to R-lambda=infinity. The models give good agreement with experiment for decaying turbulence and produce negligible SGS dissipation for forced turbulence in the limit of fully resolved flow

Geometric study of Lagrangian and Eulerian structures in turbulent channel flow

We report the detailed multi-scale and multi-directional geometric study of both evolving Lagrangian and instantaneous Eulerian structures in turbulent channel flow at low and moderate Reynolds numbers. The Lagrangian structures (material surfaces) are obtained by tracking the Lagrangian scalar field, and Eulerian structures are extracted from the swirling strength field at a time instant. The multi-scale and multi-directional geometric analysis, based on the mirror-extended curvelet transform, is developed to quantify the geometry, including the averaged inclination and sweep angles, of both structures at up to eight scales ranging from the half-height δ of the channel to several viscous length scales δ_ν. Here, the inclination angle is on the plane of the streamwise and wall-normal directions, and the sweep angle is on the plane of streamwise and spanwise directions. The results show that coherent quasi-streamwise structures in the near-wall region are composed of inclined objects with averaged inclination angle 35°–45°, averaged sweep angle 30°–40° and characteristic scale 20δ_ν, and 'curved legs' with averaged inclination angle 20°–30°, averaged sweep angle 15°–30° and length scale 5δ_ν–10δ_ν. The temporal evolution of Lagrangian structures shows increasing inclination and sweep angles with time, which may correspond to the lifting process of near-wall quasi-streamwise vortices. The large-scale structures that appear to be composed of a number of individual small-scale objects are detected using cross-correlations between Eulerian structures with large and small scales. These packets are located at the near-wall region with the typical height 0.25δ and may extend over 10δ in the streamwise direction in moderate-Reynolds-number, long channel flows. In addition, the effects of the Reynolds number and comparisons between Lagrangian and Eulerian structures are discussed